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THE EFFECT OF WALL INTERFERENCE UPON THE AERODYNAMIC CHARACTERISTICS OF AN 
AIRFOIL SPANNING A CLOSED-THROAT CIRCULAR WIND TUNNEL 

By Walter G. Vincenti and Donald J. Graham 


SUMMARY 

The results of a theoretical and experimental investigation of 
wall interference for an airfoil spanning a closed-throat circular 
wind tunnel are presented. Analytical equations are derived 
which relate the characteristics of an airfoil in the tunnel at 
subsonic speeds with the characteristics in free air. The 
analysis takes into consideration the effect of fluid compressi- 
bility and is based upon the assumption that the chord of the 
airfoil is small as compared with the diameter of the tunnel. 
The development is restricted to an untwisted, constant-chord 
airfoil spanning the middle of the tunnel. Brief theoretical 
consideration is also given to the problem of choking at high 
speeds. Results are then presented of tests to determine t he low- 
speed characteristics of an NACA 1412 airfoil for two chord- 
diameter ratios. While, on the basis of these experiments, no 
appraisal is possible of the accuracy of the corrections at high 
speeds, the data indicate that at low Mach numbers the analytical 
results are valid, even for relatively large values of the chord- 
diameter ratio. 

INTRODUCTION 

The design of modem high-performance airplanes requires, 
insofar as possible, an accurate knowledge of airfoil profile 
data at Reynolds and Mach numbers attained in flight. 
Since the size and power of wind tunnels are subject to vari- 
ous practical limitations, most existing tunnels, even if they 
can provide the desired Mach number, are not capable of 
attaining full-scale Reynolds numbers for all flight conditions. 
To minimize this shortcoming in tunnel tests of airfoil 
profiles, it is therefore necessary to use models having as 
large a chord as possible relative to the cross-sectional 
dimensions of the tunnel test section. In order to eliminate 
the effects of supporting struts and to exclude the indetermi- 
nate tunnel-boundary interference involved in the testing 
of large-chord airfoils of limited span, it has become common 
practice in such tests to use airfoils which completely span 
the test section. Even for these so-called “through” 
models, however, the tunnel-boundary interference can still 
be considerable, and accurate correction must be made for 
its effects if the tunnel data are to be used with confidence 
in the calculation of free-flight airplane characteristics. 

The t unn el-boundary interference for airfoils spanning 
wind tunnels of various types has been the subject of numer- 
ous theoretical and experimental investigations. The inter- 
ference for rectangular tunnels having rigid walls normal to 
the span of the airfoil and either rigid walls or free boundaries 


parallel to the span has been discussed theoretically by sev- 
eral writers. For example, Lock (reference 1), Glauert (ref- 
erence 2), and Goldstein (reference 3), give the necessary 
tunnel-wall corrections for an airfoil spanning a rectangular 
tunnel in an incompressible fluid; while Goldstein and Young 
(reference 4) show how these corrections, as well as those for 
any general case of interference in an incompressible fluid, 
can be modified to take account of fluid compressibility. 
Reference 5 gives the corrections for the compressible case 
in a closed-throat rectangular tunnel, as well as a critical 
discussion of the results of the previous references and some 
experimental data from low-speed tests. Fage (reference 6) 
also presents experimental drag data for several symmetri cal 
bodies of various sizes in a closed-throat rectangular tunnel. 
Experimental and theoretical results for an airfoil spanning 
a completely open-throat rectangular tunnel are given by 
Stuper (references 7 and 8). The case of an airfoil spanning 
an open-throat circular tunnel has been the subject of a 
number of investigations, including theoretical treatments 
by Glauert (reference 9), Stuper (references 7 and 8), and 
Squire (reference 10), and experimental measurements by 
Stttper (references 7 and 8) and Adamson (reference 11). 
Apparently, the case of the closed-throat circular tunnel has 
received no attention. 

Since this last case is often encountered in practice, an in- 
vestigation was made of the tunnel-wall interference at sub- 
sonic speeds for a wing spanning a closed-throat circular 
tunnel. The present paper presents the results of this in- 
vestigation. In the first part of the paper, analytical equa- 
tions are derived relating the characteristics of the airfoil 
in the tunnel with those in free air for a compressible fluid. 
Some consideration is also given to the phenomenon of chok- 
ing which occurs at high speeds. In the second part, the 
validity of the theoretical results is examined by the analysis 
of experimental data for an NACA 4412 airfoil for two ratios 
of airfoil chord to t unn el diameter. The investigation is 
restricted to untwisted constant-chord airfoils spanning the 
middle of the tunnel. 

THEORY 

As in reference 5, the theoretical development of the 
tunnel-wall corrections is divided conveniently into two gen- 
eral sections. First, the influence of the wall upon the field 
of flow at the airfoil in the tunnel is determined. Second, 
the aerodynamic characteristics of the airfoil in this field of 
flow are related to the corresponding quantities in fr ee air. 
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In this manner, simple formulas are finally obtained which 
enable the prediction of the free-air characteristics when the 
characteristics in the tunnel are known. 1 

Again, as in reference 5 , the analysis is based upon the 
method of superposition. To this end, it is assumed that the 
airfoil is of small thickness and camber, so that the. induced 
velocity is everywhere small as compared with the velocity 
of the undisturbed stream. With this assumption, the total 
induced velocity at any point is the simple vector sum of the 
separate velocities induced at that point by the interference 
between the tunnel wall and the airfoil camber, thickness, 
and wake. Thus the effects of camber, thickness, and wake 
may each be analyzed separately and superposed to obtain 
the desired result for the complete airfoil. As pointed out 
in reference 5, this procedure is permissible even in the com- 
pressible fluid if the airfoil is of small thickness and camber 
as assumed. 

Before proceeding to the actual development of the theory, 
it is useful to contrast the present problem with the problems 
of through airfoils in the various types of rectangular tunnels 
and in the open-throat circular tunnel. In the case of an 
airfoil spanning a rectangular tunnel haying rigid walls 
normal to the span of the airfoil, the problem is relatively 
simple. If the effect of the boundary layer along the tunnel 
walls is neglected, the flow is sensibly the same in all planes 
normal to the span; that is, there is clearly no spanwise 
variation in lift. The air flow is thus essentially two-dimen- 
sional, and the interference problems of camber, thickness, 
and wake can be analyzed by the customary means of a 
system of images with axe9 parallel to the span of the airfoil 
(references 1 , 2, 3, and 5). This is true whether the tunnel 
boundaries parallel to the span are fixed or free. In this 
manner, tunnel-boundary corrections can be derived for 
airfoils of moderately large chord as compared with the 
height of the tunnel test section. 

In the case of an airfoil spanning a completely free jet, 
whether rectangular or circular in section, the lift necessarily 
falls to zero at the boundary of the jet. There thus exists 
in this case a pronounced spanwise variation in lift and an 
attendant system of trailing vortices. In the existing treat- 
ments of the problem, only the interference between these 
trailing vortices and the jet boundaries is considered, the 
interference effects associated with the chordwise distribu- 
tion of bound vortices and with the airfoil thickness and 
wake being completely neglected. This procedure implies 
the assumption that the chord of the airfoil is very small 
relative to the dimensions oi the jet. In this manner, the 
problem is reduced to a limiting case of the usual problem 
of an airfoil partially spanning the jet, and, as in this latter 
case, the component of dowmvash induced at the airfoil by 
the interference between the walls and the trailing vortices 
is one-half as great as the corresponding component an 
infinite distance downstream. The theoretical determina- 
tion of the wall interference may thus be treated as a prob- 
lem of two-dimensional flow in a plane normal to the axis of 
the tunnel infinitely far behind the airfoil. The boundary 
conditions for either the rectangular or circular jet are then 
readily satisfied by the introduction of a suitable system of 
image vortices with axes parallel to the axis of the tunnel 


(references 7, 8, 9, and 10). This method of analysis, how- 
ever, is inadequate if the chord of the airfoil is even moder- 
ately large as compared with the dimensions of the jet. 

The case of the airfoil spanning a closed-throat circular 
tunnel is more complex than either of the foregoing problems. 
Unlike the condition prevailing in the free jet, the lift in 
this case need not fall, to zero at the boundary — that is, at 
the tunnel wall — so that the spanwise variation in lift is 
not necessarily large. In fact, as will be seen, the lift is 
constant across the span of the airfoil, and no system of 
trailing vortices exists. The assumption of a very small 
chord and the consequent reduction of the problem to a case 
of two-dimensional flow in a plane infinitely far downstream 
is thus without meaning. On the other hand, an analysis 
for airfoils of moderately large chord in the manner em- 
ployed in the case of the rectangular tunnel with rigid side 
walls is not possible. In the closed-throat circular tunnel 
the fiow r in all planes normal to the span of the airfoil is not 
the same, so that the effect of the bound vortices, and of the 
airfoil thickness and wake as well, cannot be treated as a 
problem in. two-dimensional flow. Furthermore, the bound- 
ary conditions at the tunnel wall cannot be satisfied for the 
actual three-dimensional problem by any known system 
of images. The solution of the problem for the closed-throat 
circular tunnel thus requires an analysis entirely different 
from those employed in the previous instances. 

INFLUENCE OF TUNNEL WALL UPON FIELD OF FLOW AT AIRFOIL 

An approach to the problem of the airfoil spanning a 
closed-throat circular tunnel is afforded by the work of von 
KArmfin and Burgers in reference 12 (pp. 206 to 273), where 
the velocity potential at an arbitrary point in the tunnel is 
determined for a U-shape vortex of infinitesimal span in an 
incompressible fluid. 



Figure I.— Elementary U -shaped vortex In closed-throat circular tunnel. 
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A system of rectangular coordinates a;, y, 2 is introduced 
as shown in figure 1. The ar-axis is taken on the center line 
of the tunnel with its positive direction downstream. The 
2 -axis is positive downward, and the y-axis positive to the 
left for an observer looking against the direction of flow. 
An alternative system of cylindrical coordinates x, a, 0 is 
defined by the relations 


y=u cos 6 
2 =« sin 6 


(1) 


The positive direction of circulation is defined so that a 
vortex with positive circulation exerts a force on the fluid in 
the direction of the positive 2 -axis. In other words, the lift 
force experienced by a positive vortex is in the negative 2 
direction. The velocity of the fluid in the undisturbed 
stream is denoted by V' and the radius of the tunnel by r. 
Other symbols are defined as introduced in the text. A list 
of the more important symbols and their definitions is given 
in Appendix C. 

Consider now a U -shape vortex of infinitisimal span d-q 
parallel to the y-axis and situated in the 2 / 2 -plane at the point 
ij=«o cos 80 , f=wo sin d a . If the strength of the vortex is 
denoted by r' the velocity potential in the closed tunnel at 
the points x, u, 6 is given by von Kfirmfin and Burgers, for 
negative values of x, as 


T’dr i f- aa 

‘ 4r 


f " to 

J» Sf^ 


( 2 ) 


where 


2 2 cos m( 6 — 8 , 

' m-0 «-l 


J ro(X,Co)<7 w (X,tOo) 

LO-^) W(X,r)J 


2(g-») 

1 s 


(3) 


(It should be noted that the quantity £ appearing in these 
equations is merely a variable of integration and has no physi- 
cal significance.) The quantity «7 m (X,«) is a Bessel function 
of the first kind of the order m. The summation with respect 
to m extends over all the positive integers and includes 
m=0; the prime added to the summation sign indicates that 
a factor % must be inserted before the term corresponding to 
m= 0. The summation with respect to s for every m 
extends over all positive roots of the equation 


Jm'M= 0 (4) 

where J m '(K,r) is the derivative of the function J m (Xp*) with 
respect to its argument. The notation used throughout this 
paper for the Bessel functions is that of Watson (reference 
13), which is the same as that of the Smithsonian tables 
(reference 14). 

By differentiating 0 with respect to f and then integrating 
with respect to |, as indicated in equation (2), the velocity 
potential becomes finally 


r,i ” ±'± 


e x ‘*J n (X,co) 


fm cos 6 0 sin m(0— 0 C ) c ^*"-* {l>ll ^ + 
L «o 

X, sin 0„ cos m(0— 0 O ) J n '(X^o 0 )J 


X 


(5) 


As pointed out, this expression applies only at negative values 
of x. As will be seen later, the necessary results for positive 
values of x can be derived from considerations of symmetry. 

By means of equation (5), it is possible to evaluate the 
wall interference associated with both airfoil camber and 
thickness for the case of the incompressible fluid. These 
results can then be modified for the effect of fluid com- 
pressibility by the methods of reference 4. It is found 
finally that, for a closed-throat circular tunnel, the effects of 
interference between the walls and the airfoil camber are 
identical with the corresponding interference effects for the 
same airfoil spanning a closed-throat rectangular tunnel, the 
height of which bears a known relation to the diameter of 
the circular tunnel. A similar conclusion is obtained regard- 
ing the effects of interference between the walls and the air- 
foil thickness, except for a numerical difference in the 
relation between the diameter of the given circular tunnel 
and the height of the equivalent rectangular tunnel. The 
interference effects associated with the wake of the airfoil 
are not analyzed in detail, but their magnitude can be 
estimated with reasonable accuracy by comparison with the 
results for the thickness effect. In order to simplify the 
complex mathematics of the problem, the interference effects 
are calculated only for the section of the airfoil at the center 
line of the tunnel. As will be seen later, however, experi- 
mental data indicate that the results are applicable at any 
spanwise station. 

Camber effect. — To analyze the effect of the interference 
between the tunnel walls and the airfoil camber, the thickness 
and wake of the airfoil are considered to be removed and the 
airfoil reduced to its mean camber line. The resulting infini- 
tesimally thin airfoil may then be replaced by a sheet of 
continuously distributed, bound vortices which, in the general 
three-dimensional case, consist of both spanwise and chord- 
wise vortices. The velocity induced at any given point on 
the camber line is then obtained by integration over the entire 
vortex sheet. As in all thin-airfoil theory the distribution 
of bound vorticity must be such that the resultant of this 
induced velocity and the free-stream velocity is tangential 
to the camber line at all points. As will be seen, however, 
the actual theoretical determination of the distribution of 
vorticity is not necessary in this case. 

In calculating the velocity field of the vortex system, it is 
assumed that the bound vorticity is distributed in the middle 
plane of the tunnel — that is, in the ay-plane — rather than 
along the camber line and that the induced velocity at any 
point on the camber line is the same as the induced velocity 
at the corresponding point in the ay-plane. From equation 
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(5), the velocity potential at any point x, u, 9 for a vortex 
element on the y-axis at the point y—v (0 o =O, wo=t?) is 


T'dy ^mainme /oX 

■~ 3 - 2 -i 2 -j > — ' W 

* m " 1 ' =1 v(l-^)x.W(X,r) 


The term for m=0 disappears by virtue of the factor m in 
the numerator of the general term. The vertical induced 
velocity v t r in the incompressible fluid is then 


/ d<t> dfl . 5<j> 5 m 
’ e = dz — d0 Dz 


For points in the ay-plane (6=0, a=y), |^=0 and ^ =— 
Thus, at points in the zy-plane, 




771-1 tf-1 


,,y(l-^)x,W(X.r) 


(7) 


The complicated double series in this equation can be 
reduced to a single series and the mathematics of the problem 
greatly simplified by limiting the discussion to the chord- 
wise section of the airfoil at the center line of the tunnel 
(y= 0, 2=0). From the known relations for Bessel functions 
(cf. reference 13), for y=0 

J\ (x«y) 1 

x«y 1/2 

^f^=0for m>l 


Thus, at points on the x-axis, 


r 'dr? eWi(X.n) 

2^ ^(l-^X^Ur) 


( 8 ) 


where the summation with respect to s extends over all the 
positive roots of the equation 

Ji'(\S)=0 (9) 

From Bessel’s differential equation 

«7i"(X^)=-(l-^p)j r 1 (X,r) .(10) 


where the double prime denotes the second derivative of the 
Bessel function with respect to its argument. Equation (8) 
can thus be written 


, T’di) ” f ,x Ji(Kv) 

’ 2rryj Xp’J 1 (X,r)J r 1 "(X,r) 


(ID 


As mentioned, this equation is valid only for negative values 
of x. 

It is apparent that the series in equation (11) is rapidly 
convergent for large negative values of x, but that the con- 
vergence is slow for small negative values. Since in the 
evaluation of the velocity induced by the vortex sheet it 


is the small values of x which are of primary importance, 
equation (11) cannot be applied directly in the present case. 
It is possible, however, by means of a method demonstrated 
by Watson (reference 15), to express the series of this equa- 
tion as combination of elementary functions and a series 
of ascending powers of x and tj. The resulting series is 
readily applicable to the present problem. 

The detailed procedure for the transformation of the series 
of equation (11) is given in Appendix A. By application 
of the final result, equation (11) may be written. 


, T'd-rj " e x ‘ z J ( X,??) T'dri r 2 -fV rx 

V ‘ ~2v ri) ^ X,rJ,(X,r) t / 1 "(X,r) = *27rn; |_ 2nj 2ijYi} 1 +i l 
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kUTc+lj 1 (2p+ 1) 
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The double summation extends over all integral values of k 
and p from zero to positive infinity. The numerical coeffi- 
cient M / 2 c*+s+n=M / a/ is given by the integral 


, _ i r- p'-m+t*) , 

(2/+l)irJo ~WW 


(13) 


Here J t (f) is a modified Bessel function of the first kind of 
order unity, and Ii(t) denotes the derivative of Ii(t) with 
respect to its argument. The numerical values of yftt for 
/=1, 2, 3, 4, are evaluated by means of a series expansion 
in Appendix A. 

It is readily shown that the first term on the right-hand 
side of equation (12) agrees with the induced velocity com- 
puted for x=0 by the more elementary theory of tunnel-wall 
interference which considers only the effects of the trailing 
vortices and their images. To this end, consider the two- 
dimensional flow in a plane normal to the axis of the tunnel 
an infinite distance downstream (fig. 2). The theory states 
that the induced velocity at a given point (y, z ) in this plane is 
twice as great as the induced velocity at the correspondingpoint 
in the plane x=0 (cf. reference 12, p. 260). In the plane 
x= oo, the trailing vortices of the U -shape vortex previously 
considered constitute a vortex pair having an infinitesimal 
spacing dt\ and situated at the point y=n, 2=0. The circu- 
lation of each vortex of the pair is V' and is directed as shown 
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Figure 2.— Section through tunnel at Infinity downstream. 
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in figure 2. The boundary condition that there shall be no 
flow normal to the wall of the tunnel can be satisfied by the 
introduction at the point 2=0, y=r*/i j of an image vortex 
pair with a spacing and with the circulation of the 

vortices directed as indicated. The vertical velocity induced 
at the tunnel center by the trailing vortex pair is 


T'dn 


vc->— is 


and the vertical velocity induced at the same point by the 
image vortex pair is 

T,/ Ml 




2x 


(O’ 


T'dn 

2xr* 


The total vertical velocity at the center of the tunnel at 
x= oo is then the sum of these two velocities; that is, 


».'(«) 


T'dn (7*+ if) 


2 tt-V 


(14) 


The vertical velocity at the center of the tunnel at *=0 is 
one-half of this value, or 


«.'(0) = 


T'dn(r*+n*) 
4 xrV 


(15) 


This value agrees with the result of equation (12) for the 
special case ®=0. Thus, the first term on the right-hand 
side of equation (12) represents the vertical induced velocity 
on the center line of the tunnel at s=0 and is attributable 
entirely to the trailing vortices and to the interference 
between these vortices and the tunnel walls. The remaining 
terms represent the variation in induced velocity due to a 
displacement a distance x upstream from the origin. These 
terms arise both from a change in the effect of the trailing 
vortices and their wall interference and from the now-active 
effect of the transverse bound vortex and its interference 
with the tunnel walls. 

Although equation (12) was deduced for negative values 
of x, it can be shown that it is applicable to positive values 
of x as well. According to von KdrmAn and Burgers (ref- 
erence 12, p. 267), the vertical induced velocity at — x is 
related to the corresponding velocity at +x by the equation 

»/ (—a) =»*'(“) — v/ (+*) 

By virtue of this relation, together with the fact that 

u/(0) =-- v/( oo ), it> follows that 

v,'(+x) -v/(0) = - [v/ (-*) -»/ (0)] (16) 

That is, the difference between the induced velocity at a 
given station x and the induced velocity at x=Q must be 
an odd function of x. The terms containing x in equation 
(12), which were derived to represent this difference for 
negative values of the variable, are seen to constitute pre- 
cisely such a function. Thus the expansion of equation (12) 
is valid for positive as well as negative values of x. 
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Figure 3. — Infinitesimally thhi airfoil spanning closed-throat circular tunnel. 


The vortex sheet which represents the entire airfoil can 
now be built up by the superposition of U -shape vortices in 
the zy-plane, and the total induced velocity found by inte- 
gration of equation (12) over the entire system. The leading 
edge of the airfoil is placed on the y-axis as shown in figure 3 ; 
the trailing edge then lies at x=c, where c is the chord of the 
airfoil. The circulation of an elementary vortex having an 
infinitesimal span dn and situated at the point £=£, y—n is 
taken to be (dT'fd%)di, where ( dT'/d £) is the vorticity per 
unit length of the chordwise section at the station y= 17. 
The vertical velocity induced at the chordwise station x on 
the center line of the tunnel by a single elementary vortex is 
given by equation (12) if x and T' are replaced by (as— £) 
and (dr'/d£)d£, respectively. The total vertical velocity 
induced by the complete airfoil is then given by the double 
integral 






r(x—£) 




(-l)V a u+, + i)^(s -{) 2 ’* 1 
£0 £0 *!(*+!)! (2p+l)I2“r 2t * +y+1) 


)J dn d£ 


(17) 


The integration of equation (17) requires a knowledge of 
(dT'/d$) as a function of n and £. Theoretically, (dr'/df) 
could be determined from the requirement that the induced 
vertical velocity at every point on the camber line must 
be such that the resultant of this velocity and the free- 
stream velocity is tangential to the camber line. This 
method of procedure leads, however, to a complicated double 
integral equation, the solution of which does not appear 
feasible. Some assumption concerning the distribution of 
vorticity must therefore be made if the problem is to be 
solved. 
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To aid in the choice of a suitable assumption, experiments 
were carried out to determine the pressure distribution, both 
chord wise and spanwise, over an airfoil spanning a closed- 
throat circular tunnel. The airfoil used in the experiments, 
which are described in detail later in the report, had an 
NACA 4412 section and was untwisted and of constant chord. 
The results of these experiemnts reveal that for such an 
arrangement the lift is sensibly uniform across the span for 
angles of attack below the stall. This fact is illustrated in 
figures 6 and 7, which show the experimental spanwise 
lift distribution for the airfoil at various angles of attack in 
wind tunnels affording chord-diameter ratios of 0.357 and 
0.625. These results were at first regarded as rather sur- 
prising. Later, however, it was realized that they are only 
what might logically be expected from general considerations 
of the conditions of flow in a closed- throat t unn el. A 
demonstration of this fact is given in Appendix B, in which 
it is shown that the lift distribution is uniform across an 
untwisted, constant-chord airfoil spanning any closed-throat 
wind tunnel, irrespective of the cross-sectional shape of the 
tunnel. Detailed examination of the pressure distributions 
from which the results of figures 6 and 7 were obtained 
reveal further that at a given angle of attack the chordwise 
pressure distribution is sensibly the same for all spanwise 
stations on the airfoil ; that is, the lift per unit chord at any 
given chordwise station is constant across the span. It is 
to be expected that this result, though obtained for a par- 
ticular airfoil, will be equally true for any ordinary camber- 
line shape. Thus it is reasonable to assume that the dis- 


tribution of bound vorticity is not a function of the spamviso 
position on the airfoil; that is. (dV'jdZ) is independent of rj. 

On the basis of this assumption, equation (17) may be 
written 



f '^/dr'\ f+ r f r 2 +7j* 

o L 


r(x-$ 


r+ 



(-l)V2u + ,4n^(g-a 2p+1 

fcl(*-H)!(2p4rf)!2“r , <* + ' +1 > 



( 18 ) 


and the integration carried out with respect to ij. The first 
two terms of the integrand, however, become infinite at the 
point 77=0. These singularities, which are duo to the effects 
of the vortices trailing from tho vortex elements on tlio 
2-axis, require that special care be taken in the integration. 
The evaluation of the integral must be carried out from — r 
to —e and from + e to -f r, and to the resulting function 
must be added the effects of the trailing vortex pairs of span 
2e which straddle the x-axis. The limit of this sum must 
then be taken as e tends to zero. The vertical velocity 
induced at the point 2 on the x-axis by the vortices trailing 
from a vortex element of span 2e symmetrically placed at 
x=£, y = 0 is 



2 (a- 

<V?+ (»-£)* 



Since theJfirst two terms of equation (18) contain only 
second-order powers of 17, the integrals from — r to — e and 
from -He to -hr will be equal. The integral of these two 
terms with respect to ij thus becomes finally 
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The integration with respect to y of the double series in 
equation (18) presents no difficulty. The expression for 
v/ thus becomes after integration 


C c (d_T/\ 

2 V I+ (V)’ , 

JoKdzJ 

L (¥) 


lent to assuming that powers of the chord-diameter ratio 
(c/d) higher than the second may be neglected in the final 
equations for the tunnel-wall corrections. The approxi- 
mation is accomplished by expanding the first term of the 
integrand in ascending powers of (2— £)/r and discarding all 
terms containing powers higher than the first and by retain- 
ing only the p=0 terms of the double series. This gives for 
the induced velocity 


•V - ' -*r-V ( l) P ^'8(t+g+I) / x Z 

^ 0 ^ 0 A:!(/;-|-l)!(2p+l)!(2/:+l)2 2 *- 1 ^ r . 
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For constant spanwise circulation, the trailing vortices 
finally disappear in the integration with respect to tj. The 
integrand of equation (19) thus represents the increment of 
vertical velocity induced by an elementary vortex of constant 
circulation completely sp annin g the tunnel. 

It will now be assumed that the chord of the airfoil is 
small enough as compared with the dimensions of the wind 
tunnel that powers of (2— £)/r greater than the first may be 
neglected in the integrand of equation (19). This is equiva- 


, 1_ C‘ d T' ) 2 

V ‘ Lrr Jo dZ ^ 

(¥) 
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dZ 


which may be written 
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By substituting the numerical values for the coeffieients- 
1) from equations (A20) of Appendix A, this equation 
may be written to an accuracy of three significant figures as 



0.579 , y r \ 
-p- (z-£) J 


( 20 ) 


The foregoing result, which was derived by assuming the 
fluid to be incompressible, can be modified for the effect of 
compressibility by the methods of Goldstein and Young. 
The modification is most readily performed by means of the 
so-called “Method II” (reference 4, pp. 5-6), which compares 
the compressible and incompressible flows for equal values of 
circulation. If the Mach number of the compressible flow 
at the position of the airfoil is denoted by M, it is readily 
shown on the basis of this method that for a given distribu- 
tion of vorticity the vertical velocity induced in a compres- 
sible fluid at any point on the center line of a tunnel of radius 
r is V 1 —AT 2 times the corresponding velocity at the same 
point in an incompressible fluid in a tunnel of radius 
r-y/l — M*. Thus, from equation (20), the vertical velocity 
v / in a compressible fluid in the actual tunnel of radius r is 


IV 


Vl-AY 4 fV^r'Xr 1 0.579 , 

J 0 (jz) ( *-£j d * (21) 


The first term of this equation represents the vertical velocity 
that would be induced by a vortex sheet of infinite span in an 
unlimited fluid field. The second term thus represents the 
interference effect of the tunnel wall. 

Equation (21) may be compared with the corresponding 
result from reference 5, which discusses the wall interference 
for an airfoil in a closed-throat two-dimensional-flow wind 
tunnel. After alteration to conform with the notation and 
sign conventions of the present paper, equation (41) of refer- 
ence 5 gives for the vertical velocity at the camber line of an 
infinitesimally thin airfoil mounted on the center line of a 
two-dimensional-flow tunnel of height h 




Comparison of equations (21) and (22) shows that an infini- 
tesimally thin airfoil spanning a closed-throat circular tunnel 
of radius r experiences at its midspan section the same inter- 
ference as would be experienced by the same airfoil in a closed- 
throat two-dimensional-flow tunnel of height 


?h= 


V6(0.579) 


r=1.686r 


or, in terms of the tunnel diameter d, 


Aj= 0.843d 


(23) 


This result makes the later deter min ation of the interference 
corrections for the circular tunnel very simple, since the cor- 
rections for the rectangular tunnel are already known. 

It is readily shown by means of equation (6) that the vortex 
system which represents the infinitesimally thin airfoil in- 
duces no axial velocity at any point in the xy-plane. It 
follows that airfoil camber has no effect upon the axial ve- 
locity or pressure gradient at the position of the model. 


Thickness effect. — The interference effects associated with 
airfoil thickness can be found by reducing the given airfoil 
to its base profile and analyzing the interference between the 
tunnel wall and this profile. The base profile is defined as 
the profile the airfoil would have if the camber were removed 
and the resulting airfoil placed at zero angle of attack. If it 
is assumed that no wake is present, the interference between 
the tunnel wall and this symmetrical airfoil can be found by 
applying the results of equation (5) to Lock’s method of 
analysis of the interference on a symmetrical body in two- 
dimensional incompressible flow. (Lock’s original analysis 
appears in reference 1; an alternative explanation of the 
method is given by Glauert in reference 2, pp. 52-57.) 

Lock's method of analysis, which assumes that the chord 
of the airfoil is small as compared with the dimensions of the 
tunnel, consists essentially in replacing the given symmetrical 
airfoil by an equivalent two-dimensional source-sink doublet 
and calculating the interference between this doublet and 
the tunnel boundaries. The strength of the doublet in any 
given case is proportioned so that it induces at a considerable 
distance from itself in free air a velocity equal to the velocity 
induced at the same point by the original airfoil. In the 
two-dimensional case, the interference flow at the position 
of the airfoil is then readily found by introducing an infinite 
series of images of the doublet such as to satisfy the condition 
that there shall be no flow normal to the tunnel boundaries 
and calculating the velocity induced at the airfoil by this 
system of images. For an airfoil spanning a closed -throat 
rectangular tunnel at mid-height, the net result of the wall 
interference for the incompressible case is to increase the 
effective axial velocity at the position of the airfoil by the 
amount 


A,F' = 


TT£ 

6A* 


where y is the strength of the doublet used to represent the 
airfoil. It is shown in references 4 and 5 that the effect of 
fluid compressibility is to increase this interference velocity 
by the factor 1/[1— (M') 2 ] s/5 , where M' is the Mach number 
of the undisturbed stream in the tunnel. Thus, in the com- 
pressible case, 


A,F'= 


za 

6A 2 [1 -(MOT 71 


(24) 


The problem of the symmetrical airfoil in a closed-throat 
circular tunnel can also be solved by replacing the airfoil by 
an equivalent doublet spanning the tunnel. In this case, 
though, the interference for the doublet cannot be found_by 
the method of images. If the doublet used is composed, 
however, of two vortices in a plane normal to the steam in- 
stead of the customary source and sink in line with the 
stream, the interference velocity can be calculated by means 
of equation (5). Since the velocity fields of the two types 
of doublets are identical, the interference calculated by mea ns 
of the vortex doublet is the same as that which would_be 
obtained if the source-sink doublet were used. 

Consider a vortex element of circulation T' and span dij 
at the point B 0 in the 3/2-plane (fig. 1). From equation (5), 
the streamwise velocity v x ’ induced at any point x, «, 6 up- 
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stream from the origin by this element and its accompanying 
trailing vortices is . 

cos 0 O sin m(0— 0 O ) 

X, sin 9 0 cos m(0— 0 O ) J n '(X,w 0 )J (25) 

At a point in the middle plane of the tunnel (ay-plane), 
6=0, a—y, and the velocity is 


»* 


T'dri 


ttT 2 


2—1 2 — 1 7 mi \ A 

m ~ 05 \ 1_ ^) x,Jm2(v) 


m cos 0 O sin m0 o ^ X ^ +X, sin 0 O cos md Q J m ' (X^> 0 ) J 

(26) 

As before, the double series in this equation reduces to a 
single series if the discussion is limited to the interference at 
the center line of the tunnel. For points on the center line, 

Jo 0^sy)—J(s{0) — 1 

Jm (X.y) = Jm(0) = 0 f or m ^ 1 

and the stream wise induced velocity becomes 


. , T drj sin 6o Jo (XfCOo) 

** = ~2^?h 


(27) 


(28) 


«/o a (X 4 r) 

From the known relations for the Bessel functions 

Jo'M = — JiO^t^o) 

so that equation (27) may be written 

_ / T'dy^ J ,z sin 0p «7, (X«mq) _ . , 9cn 

~2rrr 2 h Jo 2 (\,r) 

As required by equation (4), the summation with respect to 
s in this equation extends over all positive roots of the 
equation 

JJO^i^) — — t7i(X J r)=TT (30) 

As the next step, consider a pair of sy mm etrically placed 
elementary vortices composed of a vortex of circulation —T' 
at the point w 0 , 0 O and a vortex of circulation +r' at the 
point « 0 , 0o- From equation (29), the streamwise velocity 
induced at a point on the center line of the tunnel by this 
vortex pair and the accompanying trailing vortices is 


.. / T'drj e**® sin 0ot7j(X,wo) 


JoO^tT) 


which may also be written 




2r / tri 0 sin 0Q drj ^ d'‘ x J 1 (X t o.' 0 ) 

2 x ^0 hi J 0 *(\,r) 


(31) 


(32) 


The expression (2r / co 0 sin 0 O ) which appears in this equation 
is the product of the vortex strength and the distance between 
the vortices. 

Now let the distance between tire vortices tend to zero 
white the vortex strength increases in such a way that the 
product (2r'wo sin 0 O ) retains a constant value y. Tho 
result in the limit is an elementary vortex doublet of strength 
M and span dri at the point w® =r\ on the y-axis. The 
streamwise velocity induced on the center line of tho tunnel 
by this elementary spanwise doublet and the accompanying 
trailing vortex doublets is then 


M drj ^ J^Jii'K.rj) 
2irr J 7j hi «/ 0 *(X,r) 


(33) 


As before, the infinite series in this equation is rapidly 
convergent for large negative values of x, but the convergence 
is slow for small negative values and is nonexistent when 
r=0. Once again, however, the series can be expressed as 
a combination of elementary functions and a power scries 
which is readily applied to the problem at hand. The 
details of the transformation are given in Appendix A. 
By means of the final result; equation (33) can be written 


.. / P dr) r r 2 Npi ( 1 ) v y-Hk4.j>+i) r i ik -J v ~] 

*® ~ 2rrr* [_2(u ! +a: *) ,fl t=i ho ET( A +1)1 (2^)! J 

(34) 

where the double summation extends over all integral values 
of k and p from zero to infinity. The coefficient vnt+,+n— 
P 2 f is given by the integral 

1 f” t Sf dt 

^-(2/+i)tJ 0 im 

The numerical values of this integral for j=l, 2, 3, 4 are 
evaluated in Appendix A. 

The induced velocity for a doublet spanning the tunnel is 
now readily found by taking the doublet strength y constant 
across the span and integrating equation (34) with respect 
to t i from — r to +r. This gives finally 


(35) 


0 


„ r r 1* v N 2 ' ( 1) Va(t>y+» a;a?> 

2^ Lr 2 hi ho * 1 (* + 1 )!( 2 y) !(2it -{- 1 ) 2 Z V* 

(36) 

In the integration across the tunnel, all the trailing vortices, 
of course, disappear. 

It is apparent from the symmetry of the problem, that tho 
streamwise velocity induced by a doublet spanning the tun- 
nel must be an even function of the variable x. Equation 
(36), which was derived for negative values of x, is seen to bo 
such a function and is thus applicable to positive values of 
the variable as well. 

The values of vj for vanishingly small values of x, that is, 
at the position of the doublet, is then found from equation 
(36) by expanding the first term in ascending powers of xjr 
and discarding all terms containing second powers and 
higher and by retaining only the p = 0 terms of the double 
series. This gives 


AfS-A-s 


ifii £±«_ 


2ttt 5 |_a? 2 ho &!(&+!)! (2^+1 )2 : 


2«] 


(37) 
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After substitution of the numerical values for the coefficients 
from equations (A33) of Appendix A, tliis equation 
becomes to an accuracy of three significant figures 


| , n . 1.356,u 

' x ~ 2x2* ' 2a-? 4 


(38) 


The first term of equation (38) is the velocity induced by a 
doublet of infinite span in a field of unlimited extent. The 
remaining term therefore represents the effect of interference 
between the doublet and the tunnel wall. Thus the net 
result of the interference between the airfoil thickness and 
the tunnel wall for the incompressible fluid is to increase the 
effective stream velocity at the position of the airfoil by the 
amount 

A i F= 1|^6m ' (39) 


In any particular case, n is again equal to the strength of the 
doublet used to represent the given airfoil. 

The result of equation (39) can be modified for the effect 
of fluid compressibility by the method of reference 4. In 
this case, the modification is most conveniently performed 
by means of Method I (reference 4, pp. 3-5) which compares 
the compressible and incompressible flows for a given airfoil 
of unaltered shape and size. By this method, it is readily 
shown that the stream wise velocity induced in the incom- 
pressible fluid at any po int on the center line of a tunnel of 
radius r is (AT)* times the corresponding velocity at 

th e same po int in an incompressible fluid in a tunnel of radius 
r-y/ 1 — (M 7 ) 2 - Here M f is, as before, the Mach number in the 
undisturbed stream. The increment in axial velocity in the 
compressible case is thus 


A l V f 


1.356a 

2Tr 2 [l-(Af')T /JI 


(40) 


Comparison of equations (24) and (40) shows that, if no 
wake is present, a symmetrical airfoil spanning a closed- 
throat circular tunnel of radius r experiences at its midspan 
section the same increase in axial velocity as would be ex- 
perienced by the same airfoil in a closed-throat two-dimen- 
sional-flow tunnel of height 


h ‘-jmm r=1 - 5SSr 

or, in terms of the tunnel diameter, 

A»=0.779d (41) 

The foregoing result greatly simplifies the determination of 
the true stream conditions at the position of the airfoil in the 
circular tunnel, since the necessary equations for the rec- 
tangular tunnel are already known. 

Consideration of the symmetry of the system formed by a 
base profile spa nnin g the middle of a circular tunnel indicates 
that the intereference between the wall and the airfoil thick- 
ness does not influence the vertical induced velocity v/ at 
any point on the airfoil. Similarly, the airfoil thickness has 
no effect upon the streamwise pressure gradient in the tunnel 
at the position of the airfoil. 


Wake effect. — It is shown in general terms in reference 5 
that the interference between the wake of a body and the 
walls of a closed-throat wind tunnel gives rise at the position 
of the body to a velocity increment and a streamwise pres- 
sure gradient which are Dot present in free air. This is true 
for aDy type of body and any shape of tunnel test section. 
The magnitude of this velocity increment and pressure 
gradient m the case of an airfoil spanning a closed-throat 
rectangular tmmel can be determined approximately by 
replacing the wake by the flow from a suitable fluid source 
and the tunnel walls by an infinite system of image sources. 
In the case of the airfoil spanning a closed-throat circular 
tunnel, this treatment is no longer possible, since no system 
of image sources is known which will satisfy the boundary 
conditions at the tunnel wall. A more complex method of 
analysis could conceivably be devised for this case; however, 
since the calculation is highly approximate even in the case 
of two-dimeDsional flow, such an analysis does Dot appear 
warranted. For present purposes it is probably sufficient to 
assume that the midspan section of the airfoil in the circular 
tuDDel experiences the same velocity increment aDd pressure^ 
gradient as a result of the wake interference as does the same 
airfoil iD a rectangular tunnel of a height defined by equa- 
tion (41). This assumption leads to the simplest expression 
for the final correction to the measured drag coefficient aDd 
should give results which are reasonably accurate. If it is 
assumed that the center of the wake lies in a horizontal plane 
containing the diameter of the tuonel, it follows from con- 
siderations of symmetry that the wake interference does not 
contribute to the vertical induced velocity v/ at the airfoil. 

It has already been indicated that the interference as- 
sociated with the camber of the airfoil has no effect upon the 
stream velocity at the model. The total increase in velocity 
for the complete airfoil in the circular tunnel is thus given by 
the sum of the increments caused by the thickness and the 
wake of the airfoil. In reference 5 it is shown that for the 
analagous case of the airfoil in the rectangular tunnel, the 
true velocity V at the position of the airfoil may finally be 
written 

^ r { 1 +RW*' + TW" 1 '! (42) 

where a and r are factors dependent upon the size of the air- 
foil relative to the tunnel, A is a factor dependent upon the 
shape of the base profile, and e d is the drag coefficient of the 
airfoil as measured in the tunnel. The first correction term 
in this equation represents the velocity increment caused by 
the airfoil thickness and is found by substituting the proper 
value for the equivalent doublet strength in equation (24). 
The second correction term represents the velocity increment 
associated with the wake of the airfoil. 

The factors <r and r in equation (42) are defined by 


-SOD’ 

(43) 

il 

(44) 
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where (cjh) is the ratio of the airfoil chord to the tunnel 
height. An analytic expression for A is given in equation (3) 
of reference 5. Values of A for a number of base profiles are 
given in table I, which is reproduced from this reference. 

If it is assumed that the height of the equivalent rectan- 
gular tunnel with regard to the wake interference is the same 
as that given by equation (41) for the thickness interference, 
the true velocity in the circular tunnel is found simply by 
substituting h 2 from equation (41) for A in the factors a and r 
of equation (42). The true velocity at the midspan section 
of an airfoil spanning a circular tunnel is thus 

l /= ^'1 1 + [i — (M') 2 ] 872 A - <r2 ^~ i"— (M'y r2Ci> } 
where the factors r s and <r 2 are defined by 

r 2 = 0.3210) (46) 

and 

at= 0.339 OfiJ (47) 


A correction to the stream velocity implies corrections also 
to the stream dynamic pressure, Reynolds number and 
Mach number. . These, corrections for an airfoil spanning a 
rectangular tunnel have heen determined in reference 5 on 
the basis of the assumption that the flow . is adiabatic 
The corresponding correctipps . for tjie circular tunnel can 
be found by replacing the factors r and a in equations (29), 
(32), and (33) of reference 5 by the factors t _ 2 and at of the 
present paper. The true dynamic pressure^ q, Reynolds 
number B, and Mach number M at the midspan section in 
the circular tunnel are thus related to the corresponding 
quantities in the undisturbed stream (denoted by primes) by 
the equations . . . 




2 -(AT) 2 

'[l-(M')T /a 


A a 2 


[2-(M') J ][l + 0.4(M') 2 ] 


i w | 

(48) 


7? ff/li-i I [1-Q-7(M0 2 ][1 + 0.4(M^) 2 ] ,) 

R ~ R ( 1+ [i_(ilf , ) s ] S/iA H l-(M') 2 8 tf ) 

(49) 

kA ! i i+Q-2(M / ) 2 wmiga T , >\ 

M—M |1+ 1 — (MO 2 iCi ) 

(50) 


Numerical values of the functions of M' which appear in 
these equations are given in table II, which is reproduced 
from reference 5. 

At low Mach numbers, the terms containing t#/ in the 
equations for the corrected stream characteristics are usually 
negligible as compared with the terms containing A<r 2 . At 
liigh Mach numbers, however, where the drag coefficient is 
very large, the terms with r&t predominate. 

RELATIONS BETWEEN AIRFOIL CHARACTERISTICS IN TUNNEL 
AND IN FREE AIR 

The characteristics of the airfoil in free air are now readily 
determined in terms of the characteristics at the midspan 


section in the tunnel. It is simply necessary to apply the 
results of the preceding sections to the relations already 
derived.in reference 5 for the airfoil spanning a rectangular 
tunnel. 

Briefly, the method of reference 5 relates the section char- 
acteristics in the tunnel at an undisturbed stream velocity 
V' to the characteristics Id an unconfined stream having a 
velocity equal to the true velocity V which exists at the 
position of the airfoil iD the tunnel. The relation is ob- 
tained on the basis of equal values of the so-called cotangent 
component of lift in the tunnel and in free air, this being 
necessary to assure that the essential character of the pres- 
sure distribution over the airfoil is the same in both cases. 
By this procedure corrections are derived which may be ap- 
plied to simultaneously measured lift, moment, and drag co- 
efficients and angle of attack in the tunnel to obtain the cor- 
responding quantities iD free air. These corrections appear 
as functions of the factors A and <r, of the product rc d ', and 
of the Mach number M' of the undisturbed stream. The 
correction to the angle of attack, which arises out of the 
interference effects associated with camber, is proportional 
to a and independent of A and re/. The correction equa- 
tions for the lift and moment coefficients contain correspond- 
ing terms proportional to a alone, together with terms which 
depend upon the thickness and the wake effects and are pro- 
portional to the products A a and rc d '. The correction to the 
drag coefficient appears as two terms, proportional to A a 
and rc/, respectively. The term proportional to A a is in 
this case composed basically of two parts, one due to the 
thickness effect and one due to the wake effect. 

The correction equations for the airfoil spanning a circular 
tunnel can be derived directly by modifying the equations of 
reference 5 in accordance with the results of the preceding 
sections. Since the terms containing a exclusive of A 
appear as a result of the camber effect, the tunnel height h 
in such termB must be replaced by 0.843d as required by 
equation (23). In the terms which depend upon the thick- 
ness and wake effects and are distinguished by the products 
A a and re/, the quantity h is replaced by 0.779d in accord- 
ance with equation (41). This involves the assumption 
already mentioned that the height of the equivalent rec- 
tangular tunnel with regard to the wake effect is the same 
as that calculated for the thickness effect. 

As in reference 5, the free-air lift, quarter-cliord -moment, 
and drag coefficients referred to the true dynamic pressure q 
are denoted by the conventional symbols. The correspond- 
ing quantities measured in the tunnel and referred to the 
apparent dynamic pressure q' are denoted by the same 
symbols with primes added. The final equations for the 
corrected aerodynamic coefficients are then 

■ ft 2 -{M y 

Cl ~ Cl l 1 i-(Af') 1 FWif 74 A 
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Ci—Ci j l AOa 

\2-{Mr][l+QA{M'Y\ ,\ 
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and the corresponding angle of attack in degrees is 

“““ + 2tVi-(A/' 7 i Cl +4Cnc/i 1 

where the factor <r t is given by 

„i=0.289 (0 

and the factors r 2 and <r 2 are as already defined in equations 
(46) and (47). Numerical values of the compressibility 
factors which appear in these equations are given in table II. 
The corrected quantities correspond to the true Reynolds 
number and true Mach number as given by equations (49) 
and (50). 

From a rigorous standpoint, the foregoing corrections 
apply only to data obtained from chordwise pressure dis- 
tributions at the midspan section of the airfoil. Actually, 
as has already been pointed out in the discussion of camber 
effect, the experimental chordwise pressure distribution at 
any given angle of attack is sensibly constant across the span. 
The corrections should therefore be applicable with sufficient 
accuracy to data obtained from pressure distributions at any 
spanwise station. 

Reference 5 also includes a method for correcting experi- 
mental chordwise pressure distributions to free-air condi- 
tions in the case of an airfoil spanning a rectangular tunnel. 
The same procedure may be applied to pressure distributions 
over an airfoil spanning a circular tunnel if the factor t 
is replaced by r 2 and the factor <r by <ri wherever it appears 
alone and by <r 2 where it appears in the product Atr. 

CHOKING AT HIGH SPEEDS 

As explained in reference 5, for tests of a model in any 
closed-throat wind tunnel, there is some value of the Mach 
number M ' of the undisturbed stream which cannot be 
exceeded irrespective of the power input to the tunnel. 
This follows from the fact that at high speeds the combina- 
tion of model and wind tunnel acts essentially as a converg- 
ing-diverging nozzle, and the flow in the tunnel exhibits the 
characteristics of the flow in such a nozzle. Thus, at some 
Alach number less than unity in the undisturbed stream, 
sonic velocity is attained at all points across a section of the 
tunnel, usually in the vicinity of the model. When this 
occurs, increased power input to the tunnel serves merely 
to extend the region of supersonic flow downstream of this 
sonic section and has no effect upon the velocity of the stream 
ahead of the model. The tunnel is then said to be “choked”, 
and the Alach number M' of the undisturbed flow ahead of 
the airfoil has its maximum attainable value. This value is 
described as the apparent cho king Alach number, the word 
“apparent” being used to differentiate this value from the 
corresponding free-air Alach number M which would be 
computed from equation (50). 

842951— CO 17 


(53) 


(54) 


(55) 


If it is assumed that the section of sonic velocity is coin- 
cident with the section of minimum area between the model 
and the tunnel walls, the apparent choking Mach number 
can be obtained from elementary considerations of uni- 
dimensional adiabatic flow, as shown in reference 5. For 
the present case of a constant-chord airfoil spanning a cir- 
cular tunnel, the apparent choking Alach number M'cx if 
finally defined for air ( 7 = 1.4) by the relation 


4 

jr 



= 1 



M' ci 

6 


(56) 


where t, is the “effective” thickness of the airfoil and d is, as 
before, the diameter of the tunnel. A graph of this relation 
is given in figure 4. As a matter of interest, the results are 
shown for the supersonic as well as the subsonic flow regime. 
The region above the curve represents an impossible state 
of flow. 

In estimating the apparent choking Alach number in any 
practical case it is necessary to replace the effective thickness 
t, by the projected thickness t„ of the airfoil normal to the 
direction of flow. As indicated in reference 5, this procedure 
leads, in the case of the subsonic wind tunnel, to an over- 
estimation of AI' C * because it neglects the possible con- 
traction of a portion of the stream aft of the airfoil as well 
as the effect of the airfoil boundary layer. 
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The importance of the boundary layer and the accom- 
panying drag with regard to tunnel choking is pointed out 
in reference 6, where the apparent choking Mach number is 
calculated for a flat plate at zero angle of attack in a two- 
dimensional-flow wind tunnel. Since the projected thickness 
for the plate is zero, the unidimensional theory would indicate 
that no choking occurs. Actually, because of the fact that 
the plate experiences drag, choking does, take place. Similar 
considerations hold, of course, for a flat plate spanning a 
circular tunnel. In this case the apparent choking Mach 
number for air (7=1.4) is given by the equation 


no way is known to combine the thickness and drag effects 
in a single calculation as should logically bo done. 

It should be noted, as pointed out in detail in reference 5, 
that the flow in a tunnel at choking does not correspond to 
any flow in free air. Furthermore, for a range of Mach 
numbers just below choking, where the flow is influenced to 
any extent by the restrictions which finally promote choking, 
any wall-interference correction is of doubtful accuracy. 
This is particularly true if the model is at an appreciable 
angle of attack so that sonic velocity is attained across the 
stream on one side of the airfoil before it is on the other. 



2.8 (M' ch ) 2 


/. r i— (^>) a t> 


(57) 


A graph of this relation is given in figure 5. The effect of 
drag on choking for supersonic as well as subsonic wind 
tunnels is shown. It can be demonstrated that the points 
on the curve correspond to a Mach number of unity in the 
flow far downstream of the model where the wake has spread 
completely to the tunnel wall. Points above the curve 
represent impossible conditions of flow. In most cases 
encountered in subsonic tunnels, the apparent choking 
Mach number determined by the thickness of the airfoil and 
defined by equation (56) is usually the lower. For very thin 
airfoils at small angles of attack, however, the value of M' o, 
given by equation (57) can have the lower value. At present 



EXPERIMENT 

The experimental investigation was initiated for two rea- 
sons: (1) to determine the Bpamvise distribution of lift over 
an airfoil spanning a closed-throat circular tunnel, and (2) 
to examine the validity of the theoretical interference cor- 
rections derived in the preceding analysis. As has been pre- 
viously mentioned, the development of the theoretical rela- 
tions requires a knowledge of the variation in lift over the 
span of the airfoil. Since no theoretical or experimental 
evidence Jregarding this matter was available, the spanwiso 
variation in lift was investigated experimentally for an 
NACA 4412 airfoil for two ratios of airfoil chord to tunnel 
diameter. The results of these tests are also directly appli- 
cable to the examination of the validity of the theoretical 
correction equations. 

The experimental work was performed in a low-turbulence, 
nonreturn-type wind tunnel with interchangeable throat sec- 
tions of 14- and 8-inch diameter. The two chord-diameter 
ratios were obtained by testing the same airfoil in each throat 
section. Since the airspeed was held constant throughout 
the tests, this arrangement permitted the Reynolds number 
and the Mach number to be duplicated simultaneously for 
the two chord-diameter ratios. Id this manner the effects 
of aDy variation in these parameters were eliminated from 
the tests. 

The NACA 4412 airfoil was used because a model of suit- 
able size was already available ideally equipped for pressure- 
distributioD tests. The model, which is described in refer- 
ence 16, was of 5-inch chord and 30-inch span. This chord, 
together with the two throat diameters, gave chord-diameter 
ratios of 0.357 and 0.625. Id the tests, the airfoil extended 
through the walls of the tunnel and was clamped in tight- 
fitting support blocks which prevented any leakage of air at 
the walls. The 54 pressure orifices located around the sur- 
face of the midspan section of the model were connected to 
a multiple-tube manometer for measurement of the pressure 
distribution over the airfoil. To secure as accurate pressure- 
distributioD data as possible, alcohol was used as the manom- 
eter fluid and the liquid heights were recorded photographi- 
cally. 

Pressure-distribution records were secured at each of eight 
angles of attack from —4° to 15° at a Reynolds number of 
approximately 450,000 and a Mach number of approximately 
0.2 with the model mounted in both the 14-inch and the 8- 
ineh diameter throats. The spanwise distribution of lift was 
determined for each angle of attack by sliding the pressure 
orifices laterally from one wall to the other and recording the 
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indicated pressure distributions at a number of spanwise 
stations. The chordwise pressure distributions were plotted 
and mechanically integrated to obtain lift and quarter-chord 
moment coefficients. No drag coefficients were obtained 
because the experimental installation did not permit balance 
measurements to be made and wake surreys were not feasible. 

By testing the airfoil in both erect and inverted attitudes 
the inclination of the air stream with respect to the tunnel 
axis was determined for each throat section. The stream 
angle was found to be +0.45° for the 14-inch throat and 0° 
for the 8-inch throat. Corrections have been applied to all 
angles of attack for the measured angularity. 

The spanwise distribution of lift coefficient uncorrected for 
tunnel-wall interference is shown for the two chord-diameter 
ratios in figures 6 and 7 in which lift coefficients at various 
angles of attack are plotted as a function of the spanwise 
location of the measurement plane. 

. Curves of lift coefficient against angle of attack for the two 
chord-diameter ratios are shown uncorrected for tunnel-wall 
interference in figure 8 (a). The results given pertain to the 
section of the airfoil at the center line of the tunnel. The 
corresponding curves corrected for wall interference by 
means of equations (51) and (54) are shown in figure 8 (b). 
In applying the corrections, the term containing rc d was 
necessarily omitted as no measurements of drag were made. 
For the values of e d to be expected in such tests, however, 



Fiotoe 6.— Spanwise lift distribution for the NAC A 4412 airfoil, c/d— ,3J7. 



Fracas 7.— Spanwise lift distribution for tbe NACA 4412 airfoil. c/i-SSL 

this term would be negligible in comparison with the remain- 
ing terms so that this omission does not affect the final results. 
For purposes of comparison, section lift characteristics as 
obtained by Pinkerton from tests of a finite-span rectangular 
airfoil in the Langley variable-density wind tunnel (reference 
17) are also shown. These data correspond to an effective 
Reynolds number of 450,000 and are thus directly comparable 
to the results of the present test. 

In figure 9 (a) curves of quarter-chord moment coefficient 
against lift coefficient are shown uncorrected for t unn el-wall 
interference for both chord-diameter ratios. The same data 
are plotted in figure 9 (b) after correction for wall interference 
by means of equations (51) and (52). Also shown for com- 
parison are the corresponding data from reference 17. 

DISCUSSION 

An examinatioD of figures 6 and 7 reveals the previously 
mentioned fact that there is do appreciable variation iD lift 
over the spaD of the airfoil at all angles of attack up to those 
closely approaching the stalling angle. This observation 
holds for both chord-diameter ratios. Id the vicinity of the 
stall a spanwise variation in lift appears which becomes pro- 
gressively more erratic as the aDgle of attack is increased. 
As might be expected, this variation becomes apparent at a 
lower angle in the case of the larger chord-diameter ratio. 
The results of figures 6 and 7 corroborate the coDdusioD of 




Section pitching moment coefficient. 
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(a) Uncorrected for tunnel-wall effects. (b) Corrected for tunnel-wall effects. 


Fracas 8— Variation of section lift coefficient with angle of attack for the NACA 4412 airfoil. 



mniH in 


Section lift coefficient, c t > Section lift coefficient, c t 

(a) Uncorrccted lor tunuel trail cffcots. (b) Corrected for tunncl-wull offoels. 

Kiocee 9.— Varlallon of section pitching moment coefficient with section lift coefficient for the NACA 4412 ulrfofl. 









WALL INTERFERENCE UPON AN AIRFOIL SPANNING A CLOSED-THROAT CIRCULAR WIND TUNNEL 


239 


Appendix B for the particular case of the airfoil spanning a 
circular tunnel. 

From figure 8 (b), it is seen that the corrected lift curves 
for the two chord -diameter ratios agree almost exactly with 
0 D 6 another except at angles near the stall. Below the 
vicinity of the stall the corrected data coincide with the 
results of reference 17 except for a constant angular displace- 
ment of approximately 0.2°. In reference 16, PinkertOD 
estimates that his values for the angles of attack may be 
too large by a constant error of approximately 0.25° because 
of a possible error in the assumed direction of the stream. 
It is thought that the angles of attack of the present experi- 
ments are accurate to withiD ±0.1°. These limits of 
accuracy are sufficient to account completely for the apparent 
angular displacement. 

In the region of the stall, the corrected lift curves for the 
two chord-diameter ratios do not mutually coincide, but 
the data for the chord-diameter ratio of 0.357 agree with 
Pinkerton’s results within 2 percent. As previously men- 
tioned, Pinkerton’s tests were made with a finite-span 
rectangular airfoil, for which the cross-span variation in 
lift is necessarily large. It is not to be expected that the 
determination of ma ximum section lift from such tests 
would be as accurate as from tests of a through model, for 
which the cross-span lift variation is small. 

It is seen from figure 9(b) that the corrected moment 
curves agree satisfactorily with each other and with the re- 
sults of reference 17. 

In summary, for angles of attack below those in the region 
of maximum lift, the results presented in figures 8 and 9 
demonstrate the validity of the theoretical lift, moment, and 
angle-of-attack corrections for low Mach numbers and chord- 
diameter ratios up to at least 0.625. For angles in the vicin- 
ity of maximum lift, the corrections are not strictly applica- 
ble up to such a large chord-diameter ratio. The results of 
the present test indicate that an accurate determination of 
maximum lift can be made with a chord-diameter ratio at 
least as high as 0.35. An evaluation of the accuracy of the 
correction equations at high Mach numbers is not possible 
on the basis of the experimental evidence available at present. 
It is to be expected, however, that the maximum permissible 
chord-diameter ratios will decrease as the Mach number in- 
creases. 

The data of the present paper enable no definite conclusions 
to be drawn regarding the validity of the drag correction. 
However, in view of the accuracy of the other corrections for 
the circular tunnel and in view of the fact that the corre- 
sponding drag correction for a two-dimensional tunnel is 
known to be accurate, it is to be expected that this correction 
will give a satisfactory evaluation of the wall interference 
upon the measured drag. 

The equations of the present paper should not be expected 
to give accurate results when applied to tests in which air 
leakage occurs at the tunnel walls. In such tests the lift at 
the walls drops markedly, so that the assumption that the 
lift is uniform across the span is no longer valid. The 
importance of avoiding such leakage, if reliable airfoil 
characteristics are to be obtained, is pointed out in reference 


5 with regard, to tests in two-dimensional tunnels. The 
same general considerations apply in the case of an airfoil 
spanning a circular tunnel. 

CONCLUSIONS 

Airfoil data obtained from tests at subsonic speeds of an 
airfoil spanning the center of a closed-throat circular wind 
tunnel can be corrected to free-air conditions by means of 
the following equations: 


V=V'{l + {1 _ ( ^ /)T2 A«r, W } (45) 
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where t 2 , <ti, and a are given by 


,,-0.321 (f) 

(46) 

„=0.289(|)' 

(55) 

<r 2 = 0.339 (£j 

(47) 


and A is a dimensionless factor the value of which depends 
upon the shape of the base profile of the airfoil. (See table 
I and equation (3) of reference 5.) The remaining symbols 
are defined in Appendix C. Numerical values of the func- 
tions of M' which appear in these equations are given in 
table II. Experimental pressure distributions can also be 
corrected by proper modification of the method of reference 5 
as indicated in the text. 

Tests of an NACA 4412 airfoil at low speed for two ratios 
of airfoil chord to tunnel diameter demonstrate the validity 
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of the foregoing equations at low Mach mimbers. At angles 
of attack below the region of maximum lift, the equations are 
applicable for chord-diameter ratios up to at least 0.625, the 
maximum ratio tested. In the regioD of maximum lift a 
chord-diameter ratio of 0.35 is known , to be permissible, and 
still higher ratios may give satisfactory results. Ad examin- 
ation of the validity of the equations at high Mach numbers 
is not possible at present, but the maximum permissible 
chord-diameter ratios may be expected to decrease as the 
Mach Dumber increases. 

The tests also indicate that at low Mach numbers the span- 
wise lift distribution on an airfoil spanning a closed-throat 
circular tunDel is essentially constant except at angles of 
attack iD the immediate vicinity of the stall. This result 
corroborates the general conclusion of Appendix B, in which 


it is demonstrated that the lift is uniform across an untwisted, 
constant-chord airfoil spanning any closed-throat wind tun- 
nel, irrespective of the cross-sectional shape of the tunnel. 

The correction equations cannot be expected to apply at 
or iD the immediate vicinity of the choking Mach number, 
which is the ma ximum Mach number attainable with a given 
combination of airfoil and tunnel test section. The choking 
Mach number can be estimated by means of equations given 
m the report. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 



APPENDIX A 


TRANSFORMATION OF SERIES OF BESSEL FUNCTIONS 

The series involving Bessel functions which appear in the 
discussions of the interference effects associated with airfoil 
camber and thickness are, as pointed out in the text, poorly 
suited for use at small values of the variable x. It will be 
shown here, by means of a method demonstrated by Watson 
(reference 15), that the series may each be expressed as a 
combination of elementary functions and a convergent power 
series. The resulting series are well adapted for use in the 
present problem. The notation used for the Bessel functions 
is that of Watson (reference 13) and of the Smithsonian 
Tables (reference 14). 

Series for camber effect. — The discussion of the interfer- 
ence effects associated with airfoil camber involves the series 


TT r 6^‘ X J | (Aji;) 

1 J" (\,r) 


(Al) 


convergent for negative values of x. The summation with 
respect to s extends over all the positive roots of the equation 

J/(X,r) = 0 (A2) 


Letting j,—\r and k——x, the series may be written 


^x=S 7 






e -},Ur 
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where the summation is taken over all the positive roots of the 
equation 


<TO.)= o 


(A4) 


jfc— A T -tn. , consider the function 
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where the quantity Fi is a Bessel function of the second kind 
of order unity. This function has a simple pole at each of 
the points w— ±j, and is one-valued and analytic at all other 
points in the complex w-plane. Its residue at the point j, 
can be shown to be 


JiUtflr) 


e }, * >T 


which is identical with the general term of the series (A3). 
By the theorem of residues, the integral of the function (A5) 
taken counterclockwise around a contour inclosing the por- 
tion of the complex plane to the right of the imaginary axis 
is then equal to 2TiW x . The integral along a large semi- 
circle on the right of the imaginary axis tends to zero when 
the radius of the semicircle tends to infinity through values 
such that the semicircle avoids the poles of the integrand. 
It is thus necessary to retain only the integral along the 
imaginary axis. The contour must, however, have an in- 


dentation to the right of the origin, since the integrand has 
a pole there with residue (^-f-^/rij. If the radius of the 
indentation is made to approach zero, TFi may finally be 
written 




r*+V* 
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From the known relations for the modified Bessel functions, 
it is readily shown that 


t7i(±?Q = ±i/i(t) 

Yi(±it) = — /i(<)± — 

r/(±*)=±fi/(o+;h:/«) 

where I t and K x are modified Bessel functions of the first 
and second kind of order unity. By writing the integral 
in equation (A6) in two parts, one along the positive and 
one along the negative imaginary axis, and replacing w in 
these integrals by +it and —it, respectively, TTi then 
becomes 

w I ~£ti+ i 

Zrrj t Jo J i \l) 

or 

Wl= ~ r2 ^f + \ Jo” K ^ iv/r) sin itK/r) dt ~ 

\ J l}(t j sin dt (AS) 



The value of the first integral in this equation is given by 
Watson in reference 15 as 

f Ki(tij/r) sin (fa/r) dt— (A9) 

Jo KT 

The second integral can be evaluated by expanding the 
product Ii(tij/r ) sin (t/c/r) in ascending powers of t and inte- 
grating term by term. The series expansion for the prod- 
uct is 

«= » tfit+P+l) „2*t 1 „2p+I 

/jfo/r) sm (tKfr) =g g (&+i)i (2p+l)I 2 2t+1 r*<*+'+w 
and the term-by-term integration gives 

J Mtv/r) sin {t^fr) dt 
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The coefficient /aa+p+i >=/*'»/ is given by 

t _i. f" Em ttf dt 
M 2/_ xJo Ji'(t) 1 M 

which may be written after integration by parts 
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This integral is a constant for any given value of j. 

Reverting to the original variable x, the expansion for Wi 
may finally be written 


Wi= 
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This agrees with the result given without derivation by 
Tani and Taima (reference 18). 

For purposes of computation the coefficient n' 2f is written 


where 


M 2/ ‘ 





(A13) 
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The quantity /3 ' 2/ can then be expressed in a form suitable 
for computation by means of a method devised by Watson 
for an analagous integral (reference 15). 

As the first step, the function 


ft 

[//(f )] 2 COS (irt/b) 
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is written as a sum of partial fractions, 6 being a positive 
constant which will be fixed later. This can be accomplished 
by considering the integral 
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u+rfw 
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around the circle |w|==S in the complex plane. The inte- 
grand of the integral (A16) has poles at the points 

w=t, w—±ij„ w=± (71+1/2)6 

where j, is a positive zero of J/(w); and e=l, 2, 3, ..., 
n=0, 1, 2, 3, . . . . The residue at the simple pole at w=t 
is the function (A15). The residues at the simple poles at 
w— ±(n+ 1/2)6 are 

(— l)’*(n+l/2) 2 / 6 2/+1 


The poles &tw=±ij, are second order poles; the residues 
there are 

(-iK?v /+3 r j , _ 

(1— 17? cosh (vjt/b) (J,±ity 
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Now, the integral (Al6) taken around the circle |w|=R 
tends to zero when R tends to infinity in such a manner 
that the circle never passes through a pole of the integrand. 
It follows from the theorem of residues that the sum of the 
residues of the integrand at all its poles is zero; thus 
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By multiplying this equation by cos {vt/b) and integrating 
from — co to + co , it can be shown with the aid of certain 
integral relations given by Watson (reference 15, p. 36) that 
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and therefore 
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The firstseries in this equation converges rapidly when b is 
large, the second when 6 is small. A reasonable comijiomise 
for purposes of calcuation is to take 6=1. 

Equation (A19) with 6=1 has been used together with 
equation (A13) to determine the first four values of the 
coefficient 1 u' 2f . The final results are 


n't= — 0.999\ 
/«=- 1.627/ 
+ 6 =-9.78 
/*=— 120.8] 


(A20) 


Comparable values of n' 2 and n\ to the same number of 
significant figures are given without derivation in reference 
18. The value of in this latter reference agrees with that 
of the present paper but p\ differs by one in the third decimal 
place. The value given in (A20) has been carefully checked 
for several values of the parameter 6 and appears to be cor- 
rect. Values of n\ and + 8 apparently have not previously 
been computed. 

Series for thickness effect. — The series which appears in 
the discussion of the interference effects associated with 
airfoil thickness is 
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convergent for negative values of x. The summation with 
respect to s extends over all the positive roots of the equation 


Ji( X*r) — 0 


(A22) 


Lettingi t =X,r and k= —x as before, the series may be mitten 
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where the summation is taken over the positive roots of 




The function 
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has a simple pole at each of the points w= ±j„. Its residue 
at each of these points can be shown to be identical with the 
general term of the series (A23). Unlike the function in the 
previous series, this function is regular at the origin. Inte- 
gration around the portion of the complex plane to the right 
of the imaginary axis then gives 

I r t = 1 f“‘ £ ^ ^/r W, Mr) Y, (w) ^ 

2j-iJ_«,i 2 Jiiw) 

(A26) 

By applying the first two of equations (A7) and combining 
the integrals along the two halves of the imaginary axis as 
before, the series becomes 

ir®=- f tKi(ti}/r) cos ( tK/r)dt — — f tl&ti/r) cos ( tic/r)dt 

tJo xjo Mw 

(A27) 

The first integral can be evaluated by differentiating rela- 
tion (A9) with respect to k. This operation gives 

J o tK-SMlr) cos CA28) 

The second integral can be evaluated as before by expand- 
ing the product cos (f/t/r) in ascending powers of t 

and integrating term by term. The series expansion for the 
product is 
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and the term-by-term integration gives 


Jo ZW 11 1 C ° S dt 




(-D^awini^ 


a i=it=iE (Jc+ 1)1 (2p)! 2 st + 1 r«+^+ 1 
where the coefficient p-if is given by 
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Reverting to the original variable x, the expansion for TT* 
may finally be written 

w .=- . . . t fAsn 


2(ij 2 -(-2^) a/2 (k+ 1)! (2p)!2» +1 r“+ J ®+ 1 

The integral (A30) has been investigated by Watson (re- 
ference 15). Its value for any given r can be computed from 
the series 
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where 6 is an arbitrary positive constant. This equation has 
been used with 6=1 to determine the first four values of ptf. 
The final results are 


Ht= 0.797 
1-200 
Ms= 7.46 
/i g =96.2 


(A33) 


The first two of these values agree to the three decimal places 
with the two numerical values computed by Watson. The 
remaining two values have not previously been computed. 


APPENDIX B 

CONSTANCY OF LIFT OVER AN AIRFOIL SPANNING A CLOSED- 
THROAT TUNNEL 

Consider an infinitesimally thin untwisted airfoil of con- 
stant chord spanning a closed-throat -wind tunnel of arbitrary 
section. Such an arrangement is shown in figure 10, which 
is a section of the tunnel as seen from downstream. It is 
assumed that the flow in the tunnel is nonviscous and that 
the airfoil therefore has no drag. 

Suppose for the time being that the lift varies in some 
manner across the span of the airfoil. Any such variation 
•will be accompanied by a system of vortices trailing from the 
airfoil and extending infinitely far downstream. If the usual 
assumption is made that the trailing vortices are parallel to 
the axis of the tunnel, the flow pattern in a plane normal to 
the axis at infinity downstream must be of the nature shown 


A 


Figdtse 10.— Assumed flow pattern In plane normal to tunnel ails at Infinity downstream. 
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in figure 10. The flow pattern, in general, consists of a num- 
ber of separate sections within each of which the flow has a 
closed, circulatory character. The line AB, which represents 
the projection of the airfoil, extends across every such section, 
and each of the sections contains the filaments of a portion of 
the system of trailing vortices. The exact character of the 
flow pattern in any particular case depends upon the span- 
wise variation in lift and upon the cross-sectional shape of the 
tunnel. 

Now, consider the flow around a streamline within any one 
of the separate sections of the flow pattern — say the stream- 
line CD in the section at the left-hand side of the tunnel in 
figure 10. This streamline, like all the streamlines, inter- 
sects the projection AB of the airfoil in two points, denoted 
as C and D in the figure. The fact that in the presence of 
the tunnel walls each streamline must intersect AB in two 
points is essential to the discussion. If it is supposed for 
purposes of discussion that the direction of flow is clockwise 
as indicated, the vertical component of velocity at C is 
upward while the corresponding component at D is down- 
ward. This direction of flow corresponds to a net circulation 
in the clockwise direction for all the trailing vortex filaments 
enclosed within the streamline. 

At the position of the airfoil the pattern of transverse 
velocities induced by the trailing vortices is geometrically 
similar to the pattern at infinity downstream, only the 
magnitude of the velocities being different. Hence, at points 
on the airfoil directly ahead of point C, the vertical velocity 
induced by the trailing vortices is upward. At points 
directly ahead of point D, the velocity is downward. Thus, 
since the airfoil is untwisted, the airfoil section corresponding 
to C operates at a larger effective angle of attack than does 
the section corresponding to D. If the airfoil is of constant 
chord as assumed, this means that the lift at section C must 
be greater than the lift at section D. 

As has been pointed out, however, the trailing vortices 
discharged between sections C and D must have a net circu- 
lation in the clockwise direction in figure 10. This means 
that the circulation of the spanwise bound vortices at sec- 
tion D must be greater than at section, C. Since the 'direc- 
tion of stream flow was taken to be toward the observer, this 
in turn means that the lift at section C must be less than 
that at section D, which is in direct contradiction to the pre- 
vious result. The original supposition that the lift varies 
across the span thus leads to two mutually contradictory 
conclusions and is therefore invalid. It follows that the 
spanwise d is tribution of lift is uniform across an untwisted, 
constant-chord airfoil spanning any closed-throat wind 
tunnel, irrespective of the cross-sectional shape of the 
tunnel. 

As mentioned at the outset, this result depends upon the 
assumption that the airfoil is infinitesimally thin and has no 
drag. It will not bo strictly true if the increase in effective 
stream velocity caused by the interference between the walls 
and the airfoil thickness and wake is not uniform across 
the span. The result also neglects any effect that the bound- 
ary layer along the walls of the tunnel may have upon the 
lift distribution. That these approximations are not serious, 
at least in the case of the circular tunnel, is indicated by the 
experimental results of figures 6 and 7. 


The foregoing reasoning is, of course, inapplicable for an 
airfoil which does not span the tunnel or for a finite-span 
airfoil in free air. In these instances, the projection of the 
airfoil does not extend across all of the seetioos into which 
the transverse flow pattern is divided, and the streamlines 
of this pattern need not intersect the projection of the airfoil 
in two points. Under these conditions a type of varying lift 
distribution can be found which does not lead to a logical 
inconsistency. 
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LIST OF IMPORTANT SYMBOLS 
airfoil chord 

■ diameter of circular tunnel 
radius of circular tunnel 
height of rectangular tunnel 

height of rectangular t unn el equivalent to circular 
tunnel with regard to camber effect 
height of rectangular tunnel equivalent to circular 
tunnel with regard to thickness effect 

1/4 (j^i chord-height factor with regard to wake 
effect in rectangular tunnel 
0.321 Q0 ; chord-diameter factor with regard to wake 
effect in circular tunnel 

Q0 ; chord-height factor with regard to camber 
and thickness effect in rectangular tunnel 
0.289 00 » chord-diameter factor with regard to 
camber effect in circular tunnel 

0.339 00 > chord-diameter factor with regard to 
thickness effect in circular tunnel 
factor depending upon shape of base profile (see 
_equation (44) and table I) 
angle of attack 
section lift coefficient 
section quarter-chord-moment coefficient 
section drag coefficient 
stream velocity 
dynamic pressure 
Mach number 
Reynolds number 
rectangular space coordinates 
cylindrical space coordinates (see equations (1)) 
circulation of single line vortex in tunnel 
circulation per unit chord length 
y and z coordinates of elementary vortex 
radial and angular coordinates of elementary vortex 
chordwise coordinate of elementary vortex; also 
variable of integration in equations (2) and (3) 
velocity potential 

x and 2 components of induced velocity 
increase in axial velocity at position of airfoil in 
tunnel 

doublet strength 
projected thickness of airfoil 
effective thickness of airfoil 
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Tl't , W 2 series of terms involving Bessel functions (see equa- 
tions (Al) and (A21) of Appendix A) 

J m ,Y m Bessel functions of first and second kind of order m 
(Watson's notation) 

I m ,K m modified Bessel functions of first and second kind of 
order m (Watson’s notation) 

X, variable of s umm ation defined by the roots of the 
equation J n '(\,r)=0 

j, \,r; root of the equation J M '(j,}= 0 

P-ipp'v numerical coefficients (see equations (35) and (13)) 

k, n,p variables of s umm ation 

t,w variables of integration 

k alternate variable defined as equal to — x 

Superscripts 

(') when pertaining to fluid properties, denotes values 
in the undisturbed stream in the tunnel; when 
pertaining to airfoil characteristics, denotes values 
in tunnel, coefficients being referred to dynamic 
pressure q'; denotes first derivative of Bessel 
function with respect to its argument 

(") denotes second derivative of Bessel function with 
respect to its argument 
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TABLE I.— VALUES OF A FOR VARIOUS BASE PROFILES 



Rnnklne 

oval 

Ellipse 

Jouknwski 

Section 

Conven- 

tional 

NACA 

Sections— 

ooxx 

NACA low-drag sections 

07-0 XX 

18-OXX 

18-OXX 

18-OXX 

27-OXX 

35-OXX 

45-OXX 

63-0 XX 

64-OXX 

65-OXX 

66-OXX 

67-OXX 

a os 


0.127 


0- Ill 

0. 126 
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a 119 

a 107 
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a i(n 
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a 119 

.09 
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